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We report on a microscopic calculation of n— 3 H and p— 3 He scattering employing the Argonne 
vis and v'g nucleon-nucleon potentials with and without additional three-nucleon force. An R- 
matrix analysis of the p— 3 He and n— 3 H scattering data is presented. Comparisons are made for the 
phase shifts and a selection of measurements in both scattering systems. Differences between our 
calculation and the 7?-matrix results or the experimental data can be attributed to only two partial 
waves ( 3 Po and 3 P2)- We find the effect of the Urbana IX and the Texas-Los Alamos three-nucleon 
forces on the phase shifts to be negligible. 



Introduction 

It is well known that realistic nucleon-nucleon (NN) 
forces cannot reproduce the 3 H and 3 He binding energies. 
Three-nucleon interactions (TNIs) are added to give the 
necessary small corrections but they still fail to reproduce 
certain properties of the three nucleon system, most no- 
tably the A y analyzing power in Nd scattering jj]. Yet 
the 30% deviation of A y can be resolved by tiny changes 
in the Nd scattering phase shifts (on the order of 0.1 de- 
grees ||, §]). Furthermore very many operators can 
contribute to a TNI and the lack of stringent conditions 
in the three-nucleon system on the structure of the TNI 
makes its application to other systems desirable. In 1| it 
was shown that although a realistic NN force can gener- 
ally reproduce the 4 He system, there remain differences, 
most notably in the analyzing powers. The intensely 
studied 4 He system Q is unfortunately very difficult to 
describe due to the many resonances and the 4 He bound 
state. Therefore we investigate the much simpler systems 
p— 3 He and n— 3 H where data exist in the energy range 
of interest. 

We organize the paper in the following way: The next 
section contains a description of the i?-matrix analysis 
of the p— 3 He and n— 3 H data. After that we discuss 
the Resonating Group Model (RGM) calculation and the 
model space used. Then we compare i?-matrix and RGM 
results of the phase shifts for various model spaces, and 
demonstrate the differences for a selection of typical ob- 
servables. Finally we will discuss the effect of TNIs on 
the scattering. 



shift to correct for the short-ranged Coulomb differences 
between 4 Li and 4 H. This shift of the R- matrix eigenener- 
gies (AE\ = —0.86 MeV) for the 4 Li system was adjusted 
by hand to reproduce approximately the n — t total cross 
section measurement of Phillips, et al. lo). 

More recently, we added to the p+ 3 Kc analysis 
analyzing-power and spin-correlation data measured at 
energies between 4 and 10 MeV by Alley and Knutson 
pOfl . These high- precision data made small changes in 
the phase shifts, but did not alter the qualitative nature 
of the solution, shown as the solid curves in Figs. 1-5. 
With these data we achieve a x 2 per degree-of-freedom 
of 1.27. 

Finally, we combined the p+ 3 He data with the n — t 
total cross sections of ref. (9| and fitted both reactions si- 
multaneously. The single energy shift used earlier did not 
give a particularly good fit to the total cross sections, so 
we allowed the low-lying eigenenergies in the 4 H system 
to adjust separately from those in the 4 Li system, while 
keeping the reduced-width amplitudes the same in both 
systems. This resulted in a x 2 P er degree-of-freedom of 
1.59. Most of the increase came from the fit to the n — t 
total cross sections, which have uncertainties on the or- 
der of 0.2%. The x 2 P er point of the fit to the p+ 3 He 
data increased only from 1.22 to 1.24. The analysis is 
based on 1447 data points having proton energies be- 
tween 1.01 and 19.7 MeV and neutron energies between 
0.06 and 20.06 MeV, and allows a maximal orbital an- 
gular momentum L max = 4. The 4 Li resonance energies 
still agree with those of which also contains a brief 
description of the i?-matrix method. 



I. E-MATRIX ANALYSIS 

The i?-matrix analysis of the T = 1 part of the A = 4 
system began many years ago with an analysis of p+ 3 He 
scattering data below 20 MeV incident proton energy . 
Only one solution described all the data included, namely 
the one with 5( 3 Pi) > i5( 1 Pi), giving a x 2 P er degree-of- 
freedom value of about 1.23. Later, this solution was 
reflected to the n+ 3 H system [||, using a simple energy 



II. RGM AND MODEL SPACE 

We use the Resonating Group Model fl|, £|, |]| to 
compute the scattering in the H and Li systems us- 
ing the Kohn-Hulthen variational principle E3j. The 
main technical problem is the evaluation of the many- 
body matrix elements in coordinate space. The restric- 
tion to a Gaussian basis for the radial dependencies of 
the wave function allows for a fast and efficient calcu- 
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lation of the individual matrix elements (TT[ [li|. How- 
ever, to use these techniques the potentials must also be 
given in terms of Gaussians. In this work we use suitably 
parametrized versions of the Argonne vis and v' 8 Jl(| 
NN potentials and the Urbana IX |L6| and Texas - Los 
Alamos @ TNIs. 

The inclusion of an additional TNI requires an order of 
magnitude more computing power than the realistic NN 
forces alone. It is therefore very fortunate that enough 
data exists at low energies for the comparatively simple 
isospin T = 1 systems 4 H and 4 Li to allow for a com- 
parison between calculation and the experimental data 
or the i?-matrix analysis thereof. 

In the 4 Li system we use a model space with three 
two-fragment channels, namely the p— 3 He, the 2 H— (pp) 
and the 2 H(5' = 0) — (pp) channels. The last two are 
an approximation to the three- and four-body breakup 
channels that cannot in practice be treated within the 
RGM. The 4 Li is treated as four clusters in the framework 
of the RGM to allow for the required internal orbital 
angular momenta of 3 He or 2 H. 

For the scattering calculation we include the S, P 
and some of the D wave contributions to the — 
+ , 1 + , 0~~, 1~ and 2~ channels. From the -R-matrix 
analysis these channels are known to give essentially the 
experimental data. The full wave function for these chan- 
nels contains over 100 different spin and orbital angular 
momentum configurations, hence it is too complicated to 
be given in detail. To give an impression of the model 
space we will describe the important structures of He as 
used in the present work. 

The dominant spin S = 1/2 configurations of 3 He in- 
clude those without angular momentum and with two 
D waves coupled to total angular momentum L = 1. 
S = 3/2 in turn occurs together with a single D wave on 
each of the Jacobi coordinates or two D waves coupled to 
L = 2. Each of these configurations uses a set of one to 
three Gaussians whose width parameters were obtained 
by a non-linear optimization using a genetic algorithm 
|L8| . In this small, 29-dimensional model space we still 
achieve -6.37 MeV binding energy an rms radius of 1.78 
fm and a D state probability of 7.7% for the 3 He using 
Avis- This must be compared to — 6.92MeV known from 
Faddeev calculations |19[ whereas we find — 6.88MeV in 
a rather large model space. 

This representation of 3 He, together with a S — 1/2 
and L = 2 configuration (excited state) and the (pn) — 
(pp) fragmentations, form the basis of our 3 He scat- 
tering calculation. Once the fragment wave functions are 
fixed the scattering problem is solved with our RGM code 
relying on the Kohn-Hulthcn variational principle ifTif: 



where denotes the reactance matrix. 

The model space described above (consisting of two 
to four physical scattering channels for each J 71 ) is by 
no means sufficient to find reasonable results. So-called 



distortion or pseudo-inelastic channels |T3|] have to be 
added to improve the description of the wave function 
within the interaction region. Accordingly, the distortion 
channels have no asymptotic part. 

For practical purposes it is obvious to reuse some of 
the already calculated matrix elements as additional dis- 
tortion channels. In that way we include all the pos- 
itive parity states of the three-nucleon subsystem with 
J3 < 5/2 + in our calculation. However, it was recently 
pointed out by A. Fonseca |2^| that states having a nega- 
tive parity in the three-nucleon fragment increase the 
n— 3 H cross section noteably. Therefore we also added 
the appropriate distortion channels in a similar complex- 
ity as in the Jg 1 " case to our calculation, thereby doubling 
the size of the model space. 



III. PARTIAL WAVE ANALYSIS 

The i?-matrix analysis of the T = 1 part of the AN 
system described above represents the currently available 
p— 3 Hc and n — t data very well, as the solid lines in figs. 
Il ||and|5| show. In addition it provides elastic scattering 
phase shifts (solid line in figs. |l]and[I|) that in turn suffice 
to describe the experimental data. Therefore we take 
these phase shifts as benchmarks against which to test 
our calculations. 

In a calculation of elastic p— 3 He scattering using the 
full model space described above (including both the j£ 
and distortion channels) we find for both the Ar- 
gonne i>i8 and Wg (not shown) interactions in general 
phase shifts very similar to those given by the i?-matrix 
analysis (see figs, [l] and ^). The S- waves are negative 
due to the underlying Pauli- forbidden states, whereas 
all the P-waves are positve. The J% part increases the 
model space and therefore reduces the repulsion in the 
S'-waves and increases the (attractive) interaction in the 
P- waves. Hence, all phase shifts have to become more 
positive. The 1 5o and 3 S\ phase shifts (see fig. ||) de- 
pend only weakly on the J% part of the model space 
(long dashed line as compared to dashed line), because 
the central terms of the ./VW-interaction can only con- 
nect to a few states of the increased model space. The 
3 Po and 3 P2 phase shifts deviate markedly from the R- 
matrix results (see fig. |l|). Without the J$ components 
(long dashed line) the major difference between our cal- 
culation and the P-matrix results was the 3 P2 matrix 
element being too small, whereas the 3 Po results were 
very close. The additional distortion channels reduce 
by half the difference for 3 P2, but also raise the 3 Po phase 
considerably beyond its P-matrix values. 

It has already been pointed out previously || that the 
calculated 3 P2 matrix element in the 4 He system is too 
small (then using the realistic Bonn NN interaction pl[). 
It was shown that the manual change of 3 P2 to its R- 
matrix value removed most of the discrepancies in the 
description of the analyzing power. The part of the 
p— 3 He model space in the present calculation acts mainly 
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FIG. 1: p- 3 He P wave scattering phase shifts. RGM calculations (Avis, dashed line, restricted to postive parity in the 
three- nucleon subsystem, long dashed, and Avis + Urbana IX, dotted) compared to the 7?-matrix analysis (solid line). 



on the 3 Pq and 3 P2 channels increasing both of them by 
several degrees (fig. ^). However, the splitting of 3 P2 — 
3 Pq remains considerably too small, see fig. |?] where the 
deviations of the calaculatedp— 3 He phase shifts from the 
P-matrix values and the 3 p2 — 3 Pq splitting are shown. 
We note in passing that distortion channels as described 
above were not included in pi. 



IV. COMPARISON WITH DATA 

The increase of the 3 p2 matrix element due to the 
components suffices to describe the total and differential 
cross sections (figures || and ||) in both the n— 3 H and 
p— 3 He scattering. Again, the improvement is due to the 



larger 3 P2 matrix element, the effect of 3 Po being too 
large on the cross sections is negligible (on the order of 
1%) because of the small statistical weight. 

There remain deviations of the differential cross sec- 
tion at small and large scattering angles as well as in the 
low energy behavior of the total cross section. In the 
first case, we attribute the differences at large scattering 
angles to be an effect of the model space lacking higher 
partial waves, and our calculation agrees with the re- 
sults obtained by restricting the P-matrix partial waves 
to those included in the RGM. Even compared to these 
results, our calculation differs in the region of the inter- 
ference between the Coulomb and the strong interaction 
at small scattering angles. This discrepancy cannot be 
attributed to a single partial wave, but originates from 
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similar to the findings in [5o| . We know from |s| that the 
analyzing power is sensitive to the 3 P2 matrix element 
and it was surprising that even with the distortion 
channels increasing the 3 Pj phase shift, we cannot im- 
prove the description of A y . The reason is that most of 
the expected improvement due to the larger 3 P2 matrix 
element is canceled by the increased 3 Po phase shift that 
decreases A y , as can be seen when we lower 3 Po to its P- 
matrix value (fig. ||, where the values of A y are displayed 
close to the maximum). Increasing 3 p2 further to the 
P-matrix value will then raise A y again, but still not to 
the experimental data. Only a small part of the remain- 
ing difference is due to higher partial waves that are not 
included in the model space (see fig. ||). The remainder 
comes from small differences in the other partial waves. 



FIG. 2: p- 3 He differential cross section at P cm = 4.1 MeV. 
The experimental data is from [^] (0) and (□), the long- 
dashed line denotes the results for the Awg+TLA interaction, 
other lines as in fig. 0. 
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FIG. 3: p- 3 He analyzing power A y at P cm =4.1 MeV. The 
experimental data is from |24j (o) and fllQ[ | (A, A), the lines 
are as in fig. El. 



tiny differences between the P-matrix phase shifts and 
the calculated ones. 

The total cross section of n— 3 H scattering at low en- 
ergies depends mainly on the 3 Si phase shift whose sta- 
tistical weight 2 J + 1 is three compared to one for 1 S'o. 
Small differences in 3 Si (fig. |^ inset) lead to the 10% de- 
viation between the different models ( Avis dashed and 
long dashed, only in fig. |5|). Note that the large effect 
of the Urbana IX TNI at low energies is a combination 
of the good reproduction of the low energy 3 Si matrix 
element and the x So phase shift being too small. 

As in the 4 He system, the real challenge lies in the po- 
larization observables, e. g. the p— 3 He proton analyzing 
power A y . Even within the large model space our cal- 
culations find only about half the measured A y (fig. |^) 



V. THREE-NUCLEON FORCE EFFECTS 

Both the Av' a and Avis NN forces yield essentially the 
same phase shifts, and we convinced ourselves that the 
Gaussian parametrization of the Bonn interaction pi] ] 
also agrees with these results. In the previous section 
we showed that we cannot reproduce the p— 3 He analyz- 
ing power A y mainly due to the too small splitting of 
3 P 2 and 3 P (fig. @). Therefore we included an addi- 
tional TNI in our calculations. For Avis we chose the 
27r-exchange model implemented as the Urbana IX force 
@ and for Av' 8 the Texas-Los Alamos (TLA, (l7|) three- 
nucleon force. The latter includes only short range oper- 
ators and is claimed to resolve the A = 3 A y problem for 
a certain choice of its operator strength (ci = 3, p7|). 
In our calculation we chose this factor to be c\ = 1.5 to 
give reasonable 3 He binding energies in the model space 
used. 

In the case of the Av\% and Urbana IX interaction, the 
additional TNI leads to the expected improvement of the 
binding energy. For the small 3 He wave function we find 
-6.88 MeV instead of the -6.37 MeV for Auig alone and 
the TNI increases the D state probability Pd somewhat 
(from 7.7% to 8.2%). In a rather large model space we 
find a binding energy of —7.72 MeV and kinetic energy 
of —50.05 MeV, close to the values given in [ p)| ]. 

Both TNIs are unfortunately very time-consuming to 
evaluate and neither of them contributes significantly to 
the phase shifts (fig. [I]). Especially the Urbana IX inter- 
action does not affect the phase shifts at all while still 
improving the binding energy of 3 He. The 3 p2 — 3 Po 
splitting is somewhat improved by the TLA force but 
still far from its P-matrix values (fig. 0). Yet this im- 
provement can already be seen in the analyzing power 
(fig. ||), where the calculation including the TLA force 
gives significantly better results. Since the Urbana force 
left the phase shifts almost unchanged it affects the an- 
alyzing power (fig. |3|) only marginally. These effects can 
be more clearly seen in fig. 0, where the analyzing power 
close to the maximum is shown for the different models. 

Since the 27r exchange operators leave the scattering 
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FIG. 5: n- 3 H total cross section. The experimental data is 
from § (A), H (□) and §|] (0), the lines are as in fig. §. 
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FIG. 6: p— 3 He proton analyzing power at E C m = 4.IMeV 
and ffO degrees for different models: Avis (A), Avis J3 
only (v), Av is + Urbana IX (a), Av' 8 + TLA (♦), Av 18 
with 3 P (+) and with both 3 P and 3 P 2 (•) adjusted, R- 
matrix (pentagon) and P-matrix restricted to L < 2 (filled 
pentagon), data point (□) from jlo| ] . 

phase shifts unchanged, we considered the V% operators 
proposed in [^7j in addition to the Urbana IX interac- 
tion. The V 3 * uses a pion-exchange between the third 
particle and the center of mass of the NN system, which 
in turn uses the full N — N interaction, but only for T = 1 



x 5 3 -Si 1 P 1 3 P 3 Pi 3 Pi Z P2 - 3 Po 

FIG. 7: Deviation of the calculated p— 3 He phase shifts from 
the 7?-matrix values in the different channels and the 3 P2— 3 Po 
splitting for different models at E cul — 4.5MeV: Avis (V J% 
only, A full model), Av is + Urbana IX (a), Av' 8 (0) and Av' 8 
+ TLA (♦). 



and S = 1. The U 3 * therefore acts predominantly in the 
P-waves. In order to study the effects of this interac- 
tion within a relatively simple operator structure, we ap- 
proximated the full i-matrix by using the central part of 
the Argonne Wig potential as the NN-ioice part in U 3 *, 
and chose its strength to leave the 3 H binding energy 
unchanged while reproducing the 3 p2 matrix element at 
E cm = 2.4 MeV. 

The full model of the U 3 * NNN interaction has only a 
small influence on the binding energy of 3 H |2^] , which 
also holds for 4 He with the U 3 * potential in our calcu- 
lation. However, these operators have a large effect on 
the P-wave phase shifts. If we choose the strength of 
U 3 * so that the 3 P2 matrix element is reproduced and 
the 3 H binding remains unchanged, the 3 Pq phase shift 
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unfortunately increases again, and therefore part of the 
improvement due to the larger 3 P 2 matrix element is can- 
celled. Nevertheless, Urbana IX and V£ together achieve 
an A y of 0.35 at E cm — 4.1 MeV, a much larger effect 
than the Urbana IX force alone achieved. 



VI. CONCLUSIONS 

We have discussed a new i?-matrix analysis of the cur- 
rently available experimental data. The phase shifts cal- 
culated in this analysis were compared to an RGM cal- 
culation of p— 3 He scattering. We showed that realistic 
NN interactions describe most of the phase shifts quite 
well but fail to reproduce the 3 P 2 and 3 Po phase shifts. 
The calculated splitting between these two channels is 
much too small, and neither the Urbana IX nor the TLA 
three-nucleon force is able to improve the splitting sig- 
nificantly. In fact, it is more important to include in 



the calculation negative parity states of the three-nucleon 
subsystem than one of these two TNIs. These findings 
show that new contributions to the NNN force acting 
on the P-waves should be considered, like an LS type 
TNI, as proposed in p9] for the N — d analyzing powers, 
or the U 3 * operators proposed in |27| . 
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